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Abstract Our goal is to study the effects of the UV radiation from the first stars, quasars 
and hypothetical Super Heavy Dark Matter (SHDM) particle decays on the for- 
mation of primordial bound objects in the Universe. We trace the evolution of 
a spherically symmetric density perturbation in the Lambda Cold Dark Matter 
and MOND model, solving the frequency-dependent radiative transfer equation, 
non-equilibrium chemistry, and one-dimensional gas hydrodynamics. We con- 
centrate on the destruction and formation processes of the H2 molecule, which 
is the main coolant in the primordial objects. 
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1. Introduction 

The existence of the first objects is a direct consequence of the growth of 
the primordial density fluctuations. At the beginning, there are linear density 
perturbations which expand with the overall Hubble flow. Subsequently, these 
perturbations can grow and form primordial clouds. Clouds with enough den- 
sity contrast decouple from this flow and start to collapse. The kinetic energy 
of the infalling gas is dissipated through shocks and the cloud becomes pres- 
sure supported. The further evolution of the cloud is determined by its ability 
to cool sufficiently fast. Clouds which could not cool fast enough will stay in 
a pressure-supported stage and will not form any stars. The existence of the 
efficient cooling mechanism is necessary to continue the collapse of the cloud, 
its subsequent fragmentation and star formation. 
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In our work we are interested in the first generation of stars which are still 
forming in the low mass clouds when first luminous objects already exist. 
These objects are made from the primordial gas so they are metals free. It is 
simply because the first stars did not have much time to produce them. These 
objects could be irradiated by the UV and X-rays radiation produced by the 
first stars, quasars and hypothetical Super Heavy Dark Matter (SHDM) parti- 
cle decays (Doroshkevich & Naselsky, 2002; Shchekinov & Vasiliev, 2004). 

In the absence of metals, the most important cooling mechanism for low- 
mass primordial clouds is so called 'H2 cooling', i.e. cooling by radiation of 
excited rotational and vibrational states of H2 molecule. The presence of initial 
mass fraction of the molecular hydrogen H2 of only 10~ 6 is enough to trigger 
the final collapse of low mass clouds. 

However, molecular hydrogen is fragile and can easily be photo-dissociated 
by photons with energies of 11.26 - 13.6 eV (Lyman and Werner bands)(Shull 
& Beckwith, 1982). Destruction of the H2 would stop collapsing of the low 
mass clouds and decrease star formation rate. On the other hand X-rays ra- 
diation can increase production of the H2 by enhancement of free electrons 
fraction. From above we see that the UV and X-rays radiation background 
alters the subsequent growth of cosmic structures. It regulates star formation 
rate so it has important implication to the re-ionization history of the Uni- 
verse (Cen, 2003). It is therefore crucial to determine quantitatively the conse- 
quences of radiation feedback on the formation of early generation objects. 

Feedback of the UV background to the collapse of the spherically symmetric 
primordial gas cloud was studied by many authors (Tajri & Umemura, 1998; 
Kitayama & Ikeuchi, 2000; Kitayama, et. al., 2000; Kitayama, et. al., 2001; 
Omukai, 2001). However their calculation was simplified. They have been 
using the so called self-shielding function (Draine & Bertoldi, 1996). Our 
approach includes solving the frequency-dependent radiative transfer equation 
in the 'exact' way. We also try to check the evolution of the collapsing cloud 
in the MOND model. 

2. Description of the code 

In the Newtonian case dynamics is governed by the following equations: 
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du p dp A 

at g z dt g 

where r is the radius of a sphere of mass M, u is the internal energy per unit 
mass, p is the pressure and g is the mass density. Here eq.(l) is the continuity 
equation, (2) and (3) give the acceleration and (4) accounts for the energy 
conservation. The last term in the eq.(4) describes cooling/heating of the gas, 
with A being the energy absorption (emission) rate per unit volume. 

A consist of two parts, that is, the radiative cooling A ra ^ and the chemical 
cooling h. c hem- The former can be written as 

. dL(m) 

Kad (m) = -q ■ (5) 

The luminosity L (m) is obtained by solving radiative transfer equation. 
The net chemical cooling rate k c h em is given by 

J^chem = Q 7^ (P) 

where e c h em is the chemical binding energy per unit mass. 
We use the equation of state of the perfect gas 

p=(j- l)gu, (7) 

where 7 = 5/3, as the primordial baryonic matter after recombination is as- 
sumed to be composed mainly of monoatomic hydrogen and helium, with the 
fraction of molecular hydrogen H 2 always less than 10~ 3 . 

In case of modified gravity (MOND) we do not include the dark component 
but it is necessary to modify equation (3). Details are shown in Stachniewicz 
& Kutschera, 2005. 

In the simulations we have used the code described in Stachniewicz & 
Kutschera, 2001, based on the codes described by Thoul & Weinberg, 1995 and 
Haiman, Thoul & Loeb, 1996. This is a standard, one-dimensional, second- 
order accurate Lagrangian finite-difference scheme. 

We start our calculations at z = 500 or the end of the radiation-dominated 
era. We use our own code to calculate the initial chemical composition and 
initial gas temperature. Initial overdensities may be calculated from the power 
spectrum which may be obtained e.g. using the CMBFAST program by Sel- 
jak & Zaldarriaga, 1996. 

We apply initial density profiles in the form of a single spherical Fourier 
mode used also by Haiman, Thoul & Loeb, 1996 

6b (r) = n bQc (l + S S -^), (8) 
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where g c is the critical density of the Universe, g c = 3H 2 /8ttG with H being 
the actual value of the Hubble parameter. As the initial velocity we use the 
Hubble velocity: 

v(r) = Hr. (9) 

In our calculations we include nine species: H, H , H + , He, He + , He ++ , 
H2, H2" and e~. Full list of relevant chemical reactions, appropriate formulae, 
reaction, photoionization and photodissociation rates is given in Stachniewicz 
& Kutschera, 2001. 

3. Radiative transfer equation in spherical symmetry 

In addition to eq. (1) - (4) we solve non-relativistic, time-independent equa- 
tion for radiation transport in spherical geometry: 

du + = {ju (r) _ (r) I} (10) 

Or r on 

where I v = I u (r, //) is the intensity of radiation of frequency u, at radius r and 
in the direction = cos 9, where 9 is the angle between the outward normal 
and photon direction. j v (r) is the emissivity and k u (r) opacity at radius r. 

Our scheme of solving this equation is similar to method from Mihalas & 
Mihalas, 1984 and York, 1980. Detailed results will be described in a paper, 
which is in preparation. 
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